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Abstract 

This paper describes a self-contained, automated methodology for flow control along with 
a validation of the methodology for the problem of boundary layer instability suppression. 
The objective of control is to match the stress vector along a portion of the boundary 
to a given vector; instability suppression is achieved by choosing the given vector to be 
that of a steady base flow, e.g., Blasius boundary layer. Control is effected through the 
injection or suction of fluid through a single orifice on the boundary. The present approach 
couples the time-dependent Navier-Stokes system with an adjoint Navier-Stokes system 
and optimality conditions from which optimal states, i.e., unsteady flow fields, and controls, 
e.g., actuators, may be determined. The results demonstrate that instability suppression 
can be achieved without any a priori knowledge of the disturbance, which is significant 
because other control techniques have required some knowledge of the flow unsteadiness 
such as frequencies, instability type, etc. 
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1. INTRODUCTION 


In the last decade, increasing attention has been devoted to the development of tech- 
niques capable of enhancing our ability to control the unsteady flow in a wide variety of 
configurations such as engine inlets and nozzles, combustors, automobiles, aircraft, and 
marine vehicles. Controlling the flow in these configurations can lead to greatly improved 
efficiency and performance, while decreasing the noise levels generally associated with the 
otherwise unattended unsteady flow. Depending on the desired result, one might wish to 
delay or accelerate transition, reduce drag or or enhance mixing. There might be a need to 
postpone flow-separation, increase lift or manipulate a turbulence field. Gad-el-Hak (1989) 
and Gad-el-Hak & Bushnell (1991) provide an excellent introduction to and overview of 
various control methodologies. 

Small improvements in system performance often lead to large payoffs. For example, 
Butter (1984) estimates that a 5 percent improvement in landing maximum lift coefficient 
(Cj(max)) can translate to a 25 percent increase in payload, while a 12 to 15 percent 
payload increase could result from only a 5 percent improvement in the takeoff Cj(max). 
As further evidence of the gains of a well-controlled system, Cousteix (1992) notes that 45 
percent of the drag for a commercial transport transonic aircraft is due to skin friction on 
the wings, fuselage, fin, etc., and that a 10 to 15 percent reduction of the total drag can 
be expected simply by laminarizing the flow over the wings and the fin. This translates 
into a reduction in fuel requirements, improved performance, and/or increased payload. 
Muirhead (1978) has shown with a wind tunnel investigation that control of flow separation 
on a tractor-trailor truck can reduce the drag by 30 to 40 percent of the baseline truck 
configuration. This translates into a savings of millions of barrels of oil per year. 

Encouraged by the potential for huge rewards with what seems to be a rather modest 
input, research into ways of achieving the above gains is attaining increasing importance. In 
many technologically important situations, the flow usually starts from a smooth laminar 
state which is inherently unstable and develops instability waves downstream. These insta- 
bility waves grow expoenentially entailing nonlinear interactions, which lead ultimately to 
full developed turbulence. This process is present in various forms in most systems. There- 
fore, one goal of a good control system is to inhibit strongly, if not eliminate altogether, 
the instabilities which lead to an undesirable transitional or turbulent state. 
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1.1 The Wave-Cancellation Concept 

The simplest form of control which might achieve this objective is the wave- 
cancellation approach, based on the premise (mostly correct) that the instability mecha- 
nisms in a low-speed transition are crucially dominated by a single two-dimensional insta- 
bility wave; therefore, cancelling this wave will preclude the nonlinear interactions leading 
to laminar-turbulent transition. The wave-cancellation method further assumes that a 
wavelike disturbance can be linearly cancelled by introducing another wave equal in am- 
plitude but opposite in phase, and thus it is mostly applicable to systems governed by linear 
or quasi-linear equations. The key is to determine the parameters of the downstream wave 
which counter (cancel) the evolution of the upstream generated wave. Although there 
is little theoretical work on this topic, there are a number of experiments and numerical 
simulations which validate this approach. Both passive and active forms of control are 
possible, and an active control might or might not make use of feedback. Passive con- 
trol makes use of gemetrical and physical characteristics to affect the flow. There is no 
wave-dependent input by an external device which might require energetic input. In the 
passive case without feedback, one inputs a wave with the proper phase and amplitude 
downstream of the source (actuator), using knowledge of the upstream input. Better and 
more robust results can be achieved via active control with feedback. This technique iden- 
tifies the characteristics of the wave downstream of the actuator using a sensor, and feeds 
back (time-dependent) information to the actuator. Different types of controls, all leading 
to different levels of wave-cancellation have been studied, both experimentally and numer- 
ically. The aim is of course to achieve the perfect cancellation, wherein the input wave has 
the same shape as the input wave at one or several locations on the body, and exactly 7 r 
radians out of phase. Of course, this is not achievable in practice, as system imperfections 
and nonlinearities prevent perfect wave cancellation. 

Most of the experiments, aimed at verifying the wave-cancellation concept, were con- 
ducted on the flat plate, except those of Ladd & Hendricks (1988), Pupator & Saric 
(1989), and Ladd (1990), who considered axisymmetric bodies. Many of these experiments 
were conducted in water tunnels. Vibrating wires (Milling 1981), hot strips (Liepmann & 
Nosenchuck 1982a, b), suction and blowing (Pupator & Saric 1989; Ladd 1990), electromag- 
netic generators (Thomas 1983), and adaptive heating (Ladd & Hendricks 1988), are some 
of the methods that were used to generate the controlling wave. All these input mecha- 
nisms give the experimenter control over the phase and amplitude of the input wave. The 
experimentalists met with varying degrees of success. Among the more successful studies 
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were those of Milling (1981) and Thomas (1983) who achieved at least an 80 percent re- 
duction in the input amplitude of the 2-D wave (with 0.6%-l% amplitude). However, it 
was not possible to achieve relaminarization, probably because of the three-dimensionality 
of the flow resulting from the interaction between background disturbances and the pri- 
mary 2-D wave. As expected, the studies conducted on axisymmetric bodies produced 
relatively less wave-cancellation since these flows are highly three-dimensional. Further- 
more, good wave-cancellation results requires a linear system with constant coefficients. 
This requirement is clearly violated for a flow over a body with curvature. 

In addition to the aforementioned experiments, several theoretical (i.e. , linear com- 
putations and theory) and computational studies (i.e, nonlinear simulations) have focused 
on understanding the physics of this wave-cancellation process. Maestrello & Ting (1984) 
provided a linear asymptotic analysis to demonstrate the relationship between the input of 
localized disturbances and their effect on the Tollmien-Schlichting (TS) instability waves 
present in the wall-bounded shear flow. They showed that small amounts of local peri- 
odic heating could excite disturbances that actively control the TS waves which propagate 
in a boundary layer on a flat plate. Analogous to the experiments, several wave input 
mechanisms were considered. In one of the early Navier-Stokes simulations of active con- 
trol, Biringen (1984) used suction and blowing at the wall in a channel flow. He observed 
approximately a 50 percent reduction in the amplitudes of the 2-D instabilities and a de- 
crease in the growth of the 3-D instabilities. The Reynolds stress originally generated 
by the waves was all but removed. On the other hand, Metcalfe, Rutland, Duncan & 
Riley (1985) studied the effect of a moving wall on unstable waves traveling in a laminar 
flow on a flat plate. The simulations were based on the Navier-Stokes equations solved 
within a temporal framework. An energy analysis revealed that the wall motion causes 
the Reynolds-stress term to become negative, which implies a feed of energy from the 
perturbed flow back into the mean flow. In effect, this energy analysis showed how a per- 
turbation to an unstable flow can be stabilizing. However, an instability wave eventually 
formed downstream of the control, with the same growth rate as the uncontrolled wave. 
This is a clear indication that the cancellation was not complete. Although intuitively 
obvious, until the work of Bower, Kegelman, Pal & Meyer (1987) and Pal, Bower & Meyer 
(1991), it was not known that perfect cancellation could be obtained within the context of 
linear theory (for which the mean flow is independent of the propagating direction. They 
used the 2-D Orr-Sommerfeld equation to study and control instability-wave growth by 
superposition, and showed, within the limits of linear stability theory and the parallel-flow 
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assumption, that single and multifrequency waves can be cancelled. The basic conclusions 
drawn by the early experimentalists were confirmed by the studies of Laurien & Kleiser 
(1989) and Krai & Fasel (1989). They showed that it was possible to delay or accelerate 
(but not eliminate) transition by appropriately (either by heating or by suction and blow- 
ing) superposing disturbances out of (in) phase with the primary TS wave. To get the best 
results required that the control be applied in the linear stages of transition, before the 
secondary instability sets in. Similar results were also reported by Danabasoglu, Biringen 
& Street (1991). 

All of the previous active-control studies were undertaken with the a priori assump- 
tion that wave cancellation was accomplished by the linear superposition (or forcing) of 
waves with 7 r radian phase shifts, or one-half wavelength/period phase shifts. None of 
these previous studies were able to achieve complete (or exact) instability removal (wave 
cancellation) from the flow, except for the linear studies reported by Milling (1981) and 
Pal, Bower & Meyer (1991), where the limitations of hydrodynamic linear instability the- 
ory were assumed. In practice, complete cancellation is not possible with a single actuator 
because any small residual of the cancelled wave will grow exponentially with downstream 
distance. Therefore, the wave will reappear some small distance downstream. 

Although there have been quite a few studies on the so-called wave-cancellation ap- 
proach, no effort was made to demonstrate clearly that wave cancellation was in fact 
responsible for the observed reduction in wave amplitude. Recently, Joslin, Erlebacher & 
Hussaini (1995) performed a numerical experiment which served to unequivocally demon- 
strate the link between linear superposition and instability suppression. To ensure that 
linear superposition of individual instabilities was, in fact, responsible for the results found 
in previous experiments and computations, they carried out three simulations with i) only 
the disturbance; ii) only the control; and iii) using both disturbance and control, which 
is the wave-cancellation case. By discretely summing the control-only and forcing-only 
numerical results, they found that this linear superposed solution is identical to the wave- 
cancellation results. These tests deary verify the hypothesis that linear superposition is 
the reason for the previous experimental and computational results. 

From the above experiments, linear computations and nonlinear simulations, several 
common features emerge: i) It is impossible to achieve perfect wave cancellation unless the 
system is linear, with constant coefficients; ii) the efficiency of wave cancellation decreases 
as the system becomes more nonlinear; iii) as the geometry of the configuration becomes 
more complex, cancellation becomes more difficult; and iv) the current approaches require 
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foreknowledge of the instability wave characteristics, such as its frequency and amplitude. 

1.2 Optimal Control Theory 

The optimal control theory provides an approach which does not require a priori 
knowledge of the flow characteristics. The goal of optimal control theory is indeed to 
minimize or maximize an objective function in a robust manner. When the flow is time- 
dependent, and a strong function of initial conditions, it becomes difficult to establish the 
precise controls that will achieve the desired effect. Wave-cancellation, as discussed above, 
only works well when the input wave has a dominant frequency, and its properties are 
known. Then (either in a passive fashion, or through a feedback mechanism), one seeks to 
cancel its effect while still in a linear regime. In practice, there are many waves, which can 
interact nonlinearly in ways not always known in advance. Rather than try to cancel the 
incoming waves, one seeks appropriate controls in other ways. One means of achieving this, 
without an extensive search over the space of possible controls, is to postulate a family 
of desired controls (e.g., an arbitrary time-dependent amplitude and a specified spatial 
distribution), and an objective function (i.e. , stress over a region of the plate). Then, 
through a formal minimization process, one derives a set of differential equations, and 
their adjoints, whose solution produces the optimal actuator profile (among the specified 
set). While the solution to this set of equations cannot be accomplished in real time, the 
results can be applied using standard passive or active control mechanisms. The advantage 
is that entire collections of controls can be studied simultaneously, rather than one at a 
time. 

Optimal control methodologies have been recently applied to a variety of problems in- 
volving drag reduction, flow and temperature matching, etc. to provide more sophisticated 
flow control strategies in engineering applications. This is possible because computational 
fluid dynamics (CFD) algorithms have reached a sufficiently high level of maturity, gener- 
ality, and efficiency so that it is now feasible to implement sophisticated flow optimization 
methods, which lead to a large number of coupled partial differential equations. Optimal 
control theory is quite mathematical, and its formal nature is amenable to the derivation of 
mathematical theorems related to existence of solutions and well-posedness of the problem. 
Two recent surveys of the mathematical theories of optimal flow control are Gunzburger 
(1995) and Borggaard et al. (1995). For a mathematical treatment of a problem very 
similar to the one considered in this paper, one may consult Fursikov, Gunzburger & Hou 
(1995). Optimal control techniques will not provide the real time control that one is ul- 
timately interested in, but by systematically computing the best control within specified 
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tolerances, with a given objective function, it will be possible to develop strategies (active 
or passive) to control a wide variety of disturbances). For example, to effectively control 
a boundary-layer transition due to the interaction of cross-flow and TS wave using pe- 
riodic heating and cooling, optimal control would allow: 1) a determination of the best 
objective function to use for a given type of control (some are better than others), and 
2) provide insight into the relationship between the time dependence of the control and 
the input waves. This insight could then be built into a neural network, or other type of 
self-learning system, to allow effective control over a wide range of input parameters. 

1.3 The Current Approach 

The methodology of the current paper is based on defining a control mechanism and 
an objective for control, and then finding, in a systematic and automated manner, controls 
that best meet the objecive. In the present setting, an objective or cost functional is 
defined that measures the “distance” between the measured stresses, and their desired 
values along a limited section of the bounding wall and over a specified length of time. 
One may interpret the objective functional as a “sensor,” i.e., the objective functional 
senses how far the flow stresses along the wall are from the corresponding desired values. 
To control the flow, we imposed time-dependent injection and suction along a small orifice 
in the bounding wall. Although the spatial dependence of the suction profile is specified (for 
simplicity), the optimal control methodology determines the time- variation of this profile. 
However, unlike feedback control methodologies wherein the sensed data determines the 
control through a specified feedback law or controller, here the time-dependence of the 
control is the natural result of the minimization of the objective functional. This scenario 
is shown in Figure 1. We have a sensor that feeds information to a controller that in turn 
feeds information to the actuator. However, in the optimal control setting, the sensor is 
actually an objective functional and the controller is a coupled system of partial differential 
equations that determine the control that does the best job of minimizing the objective 
functional. The present active-control approach is demonstrated for the evolution and 
automated control of of spatially growing 2D disturbances in a flat-plate boundary layer. 
As the length of time over which the minimization process is increased, we recover the 
results obtained by wave cancellation, thus validating the approach. The ultimate goal 
of this line of research is to introduce automated control to external flows over realistic 
configurations such as wings and fuselages, and to devise novel flow control techniques. 
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2. THE OPTIMIZATION PROBLEM 


2.1. The state equations 

Let Q denote the flow domain which is the semi-infinite channel or boundary layer 
\x > 0, 0 < y < h\, where h is the location of the upper wall for the channel or the 
truncated freestream distance for the boundary layer. Let T denote its boundary and let 
(0, T) be the time interval of interest. The inflow part of the boundary \x = 0, 0 < y < h] 
is denoted by T^ and the part of the boundary on which control is applied (i.e., along 
which the suction and blowing actuator is placed) by T a which is assumed to be a finite 
connected part of the lower boundary (or wall) \x > 0, y = 0]. Solid walls are denoted 
by r w ; for the channel flow, T^, is the lower boundary \x > 0, y = 0] with T a excluded 
and the upper boundary \x > 0, y = h\; for the boundary layer flow, T^, is only the lower 
boundary with T a excluded. For the boundary-layer case, the upper boundary \x > 0, 
y = h], which is not part of T^,, is denoted by r e . Controls are only activated over the 
given time interval t £ (To, Ti), where 0 < To < T\ < T. 

The flow field is described by the velocity vector (u, v) and the scalar pressure p and 
is obtained by solving the following momentum and mass conservation equations 


du du du dp d f du\ d f du dv\ 
dt + dx + dy + dx dx\dx) dy\dy~^dx) 
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in (0, T ) , 
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Here, the initial velocity vector (uo(x,y),vo(x,y)) and the inflow velocity vector 
(ui(t,y),Vi(t,y)) are assumed given and the base flow is assumed to be Poiseuille flow 
for the channel case and Blasius flow for the boundary-layer case. The above system holds 
for both the channel and Blasius flow cases; in the latter case, the upper boundary is not 
part of r w and the additional boundary condition 


u\r ( 


= U a 


and 


p — 2v 


dv 

dy 


= P n 


m 


(0 ,T) 


( 9 ) 


is imposed, where U 0 Q and P ^ denote the free-stream flow speed and pressure, respectively. 

The particular form of the viscous terms in (1) and (2) is necessary in order to correctly 
interpret some boundary integrals that will appear below. 

The control functions gi(t, x) and < 72 ^, x) which give the rate at which fluid is injected 
or sucked tangentially and perpendicularly, respectively, through T a are to be determined 
as part of the optimization process. In order to make sure that the control remains bounded 
at To, it is required that 


gi\t=T 0 = gio(x) and g 2 \ t=To = g 20 (x) on T a , 


( 10 ) 


where g 10 (®) and g 2 o(x) are specified functions defined on T a . Commonly, one chooses 
gw(x) = g 20 (x) = 0. 

2.2. The objective functional and the optimization problem 

Assume that T s is a finite, connected part of the lower boundary \x > 0, y = 0] which 
is disjoint from T a and that (T a ,T&) is a time interval such that 0 < T a < T& < T. Then, 
consider the functional 
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where g\ and g 2 denote the controls and T\(t, x) and T 2 (t, x) are given functions defined on 
(T a , Tb) X T s . Note that since T s is part of the lower boundary of the channel or boundary- 
layer wall, vdu/dy and — p + 2 vdv/dy are the shear and normal stresses, respectively, 
exerted by the fluid on the bounding wall along T s and thus t\ and T 2 may be interpreted 
as given shear and normal stresses, repectively. Then, the boundary segment T s can be 
thought of as a sensor which measures the stresses on the wall. Thus, in (11), T s is the 
part of the boundary T along which one wishes to match the shear and normal stresses to 
the given functions t\ and T 2 , respectively, and (T a , T&) is the time interval over which this 
matching is to take place. (There are no diffllculties, other than notational, introduced if 
one wishes to match each component of the stress vector over a different boundary segment 
and/or over a different time interval.) 

The third and fourth terms in (11) are used to limit the size of the control. Indeed, no 
bounds are a priori placed on g\ or </ 2 ; their magnitudes are limited by adding a penalty to 
the stress matching functional defined by the first two terms in (11). The particular form 
that these penalty terms take, i.e, the third and fourth terms in (11), is motivated by the 
necessity to limit not only the size of the controls g\ and g 2 , but also to limit oscillations. 
The constants op, « 2 , /?i , and f3 2 can be used to adjust the relative importance of the 
terms appearing in the functional (11). 

The (constrained) optimization problem is given as follows: 


find u,v ,p, g\, and g 2 such that the functional J(u,v , p, g\, g 2 ) given 
in (11) is minimized subject to the requirement that (l)-(8) and (10) are 
satisfied and, for the boundary -layer flow case, (9) is also satified. 


3. THE OPTIMALITY SYSTEM 

We first consider, in Sections 3. 1-3.5, the case of a channel flow; the optimality system 
for the boundary layer flow is considered in Section 3.6. 


3.1. The Lagrangian Functional 

The method of Lagrange multipliers is formally used to enforce the constraints (l)-(3) 
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and (5). To this end, the Lagrangian functional 
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( 12 ) 


is introduced. In (12), u and v are Lagrange multipliers that are used to enforce the x 
and y-components of the momentum equation (1) and (2), respectively, p is a Lagrange 
multiplier that is used to enforce the continuity equation (3), and si and s 2 are Lagrange 
multipliers that are used to enforce the x and y-components of the boundary condition 
(5), respectively. Note that Lagrange multipliers have not been introduced to enforce the 
constraints (4), (6)-(8), and (10), so that these conditions must be required of all candidate 
functions u, v, p, g i, and y 2 . 

Through the introduction of Lagrange multipliers, the constrained optimization prob- 
lem is converted into the unconstrained problem: 

find u, v, p, gi, y 2; u, v, p, si, and s 2 satisfying (4), (6)-(8), and (10) 
such that the Lagrangian functional £(u, u, p, yi , y 2 , h, D, p, si , s 2 ) 
given by (12) is rendered stationary. 

In this problem, each argument of the Lagrangian functional is considered to be an inde- 
pendent variable (only subject to the constraints (4), (6)-(8), and (10)) so that each may 
be varied independently. 

The first-order necessary condition that stationary points must satisfy is that the first 
variation of the Lagrangian with respect to each of its arguments vanishes at those points. 
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One easily sees that the vanishing of the first variations with respect to the Lagrange 
multipliers recovers the constraint equations (l)-(3) and (5). Specifically, 
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where 8C/8u denotes the first variation of C with respect to it, etc. 

3.2. The Adjoint Equations 

Next, set the first variations of the Lagrangian with respect to the state variables it, 
i>, and p equal to zero. These result in the adjoint or co-state equations. Note that, since 
for the channel flow, candidate solutions must satisfy (4), (6)-(8), and (10), one has that 
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where n\ and n 2 denote the x and y components, respectively, of the outward normal to 
0 along r. Choosing variations Sp that vanish on the boundary T but which are arbitrary 
in the interior 0 of the flow domain yields that 


du dv 
dx dy 


0 on (0, T) x ft. 


( 14 ) 


Now choosing variations Sp that are arbitrary along the boundary T yields that 
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where T\r s denotes the boundary T with deleted. We note that in the above derivation 
of (14) and (15), as in the derivations found below, the boundary integrals at infinity do 
not make any contribution due to the fourth equation of (13). 
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Applying Gauss’ theorem enough times to remove all derivatives from the variation 8v in 
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the integrals on O yields 
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r- T 


Sv 


o Jr 


a 


du dv dv 

s 2 + pn 2 + vuni + vvn 2 + v— n 1 + v— n 1 + 2 v—n 2 

dy dx dy 


dYdt 


— v 


■ T 


dSv 

v — — ni + zi 


dSv 


„ \ 


0, 


where we have used (13) to eliminate boundary integrals along I\, r w and as x — > oo and 
an integral over O at t = 0. First, variations Sv that vanish at t = 0, t = T, and in a 
neighborhood of T are chosen, but which are otherwise arbitrary. Such a choice implies 
that all boundary integrals in (16) vanish giving 


dv ^du ^ dv dv dv dp 

dt + dy ^ dy dx dy dy 

d ( du dv\ d ( dv\ 
dx \dy + dx ) dy\dy) 


in (0, T) x (1 


(17) 


where equation (3) is used to effect a simplification. Next, variations that vanish in a 
neighborhood of T, but which are otherwise arbitrary, are chosen to obtain 


v\ t =T = 0 in 11. 


(18) 


Now, along T, Sv and dSv/dn may be independently selected, provided that (13) is satisfied, 
where d/dn denotes the derivative in the direction of the outward normal to O along T. If 
Sv = 0 and dSv/dn varies arbitrarily along T, then 


0 on < 


v = < 


f (o, t) x r\r, 
(o, T a ) x r s 
((T b ,T)xT s 

I r, d V 

a 2 [ -p + 2 ;/ — ~ r 2 


on (T a ,T b ) x r. 


(19) 
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To see this, note that along the inflow, I\, n 2 = 0 and d/dn = —d/dx while along the top 
and bottom boundaries n\ = 0, d/dn = ±<9/<9y, respectively, and, since 8v = 0, dSv/dx = 
0. Note that (15) and (19) agree on the boundary segments where they simultaneously 
apply. Finally, 8v is arbitrarily chosen along T a to obtain 


S2 


( du dv \ dv 

= -pn 2 - v(uni + vn 2 ) - ui— + — jni -2u—n 2 on (0,T)xT a . (20) 


Next, consider 8C/8u = 0 which yields 

f Tb f dSu ( du 
-Oil / / t— — V- T! ) dT 


T a Jr, dy V dy 

f T f f dSu dSu du dSu d 2 Su d 2 Su 

+ u \ T + 1- 8u — + V— 2v-^-z V „ „ I dildt 


J o jq 

+ [' f ( rra- 
J o Jn 


dt dx dx dy dx 2 dy 2 

dv „ d 2 8u „dbu\ f T f 

— w — — 1- p — — dildt + / / s\8u dL dt = 0 . 

dx dxdy dx J 


>o Jr„ 


Applying Gauss’ theorem enough times to remove all derivatives from the variation 8u in 
the integrals on O yields 

f Tb f dSu f du \ 

-ai / / t— — v- ti) dr 


T a Jr, dy V dy 

f T f ( du „ du ^dv „ du ^dv 

+ / / ° u u— + u— + v — 

Jq J £1 ^ dx dy dx dx 

du du dp d 2 u d 2 u d 2 v \ 

-u- v- 2 v—— - v — — - v 1 dildt 

dx dy dx dx 1 dy z dxdy 

+ f (u6u)\t = TdQ 


in 


r-T 


/ „ „ „ du du dv\ 

+ / / bu si + pni + uuni + uvn 2 + 2u——ni + v--n 2 + v——n 2 didt 

Jo Jr a \ dx dy dx J 

[ T [ ( „ dbu ^ dbu ^ dbu \ 

— v / / 2 u— — ni + u— — n 2 + v— — n i ) didi = 0 . 

Jo Jr \ dx dy dy J 

Applying to this equation the same process that led from (16) to (17)-(20) yields 


du ^du ^ dv du du dp 

dt dx dx dx dy dx 

d f^du\ d f du dv\ . 

~ U dJ\ 2 dj) ~ U d^\d^ + dJj = ° m (°’ ' 


( 21 ) 
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u = < 


and 


si = 


u\t=T = 0 in O, 

' (0,T) x T\T S 

0 on < (0, T a ) xT s 
AT b ,T)xT s 
_ Tl ) on (T a ,Tb) x T s , 


Ou ( Ou Ov \ 

-pnx - u^unx + vn 2 ) - 2V—UX - vl— + — \n 2 on (0,T)xT o 


( 22 ) 


(23) 


(24) 


In deriving (23) we have used the assumption that T s is part of the lower boundary of the 
channel so that along T s we have that ri2 = — 1. Again, there is no conflict between (15) 
and (23) along boundary segments on which both apply. 

3.3. The Optimality Conditions 

The only first-order necessary conditions left to consider are 6C/8g i = 0 and 8C/8g 2 = 
0. (These conditions are usually called the optimality conditions.) Now, since all candidate 
functions g\ and g2 must statisfy (10), it follows that 8g\ = 0 and 8g2 = 0 at t = To. Then, 
8C/8g 2 = 0 yields that 


-T, 


fh 


'To 


dg2 d8g 2 
dt dt 


<-T 


+ ^28 g 2 )dVdt + 


S2 8g2 dTdt = 0 


'To 


Applying Gauss’ theorem to remove all derivatives from the variation 8g2 yields 



d 2 g 2 1 

0,9 + 92 + ~ 7 T S 2 

dt 1 p 2 


dTdt + J (tg2-j^- 


dT 

t=Ti 


0, 


where the fact that 8g2\t=T 0 = 0 has been used. Choosing variations 8g2 that vanish at 
t = T\ but which are otherwise arbitrary yields 

d 2 g 2 1 

— 8 u 2~ 92 = —~rS2 on (T 0 ,T 1 )xT a 

at z P2 

or, using (20), 

+n = ~h{ i,+2v %) 011 (T »' Ti|xr “’ (25 > 
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where (19) and the assumption that T a is part of the lower boundary so that, along T a , 
n\ = 0 and ri 2 = — 1 have been used. Now, choosing variations that are arbitrary at t = Tj 
yields that dg 2 /dt = 0 along T a at t = Tj so that, invoking (10), g 2 (t,x) satisfies 

= 0 on T a . (26) 

t=Ti 

Note that, given p and 7, (25)-(26) constitute, at each point x £ T a , a two-point boundary 
value problem in time over the interval (To,Ti). 

In a similar manner, setting 6C/8g i = 0 yields that 


92\t=T 0 = 92o(x) and 


092 

dt 


d 2 g 1 


1 / du 

P i \ dy 


on (T 0 ,Ti) x T a 


(27) 


and 

9i\t=T 0 = 9io(x) and 

3.4. Finite Computational Domains 

In the computations, the semi-infinte domain O (we are still only considering the 
channel flow case) is replaced by a finite domain Qc defined by the introduction of the 
outflow boundary T 0 given by \x = T, 0 < y < h\. Thus, we have that Qc is the rectangle 
[0 < x < T, 0 < y < h\. We treat the ouflow by replacing (1) and (2) with 


dgi 


dt 


= 0 


on T n 


t=T, 


(28) 


du du du dp , d 2 u d 2 u „ . ^ 

W + + x = 0 m (0 5 r )x^c, (29) 

dt dx dy dx dx z dy z 


and 


dv dv dv dp d 2 v 

&+“&+W + ar !/a(l) a? 


d 2 \ 


— v 


dy 2 


= 0 in (0, T) x Oc, 


(30) 


respectively, where a(x) is a smooth function that is unity in the bulk of the flow, that 
vanishes in a neighborhood of the outflow boundary T 0 , and which is smooth throughout 
the flow. This treatment of the outflow does not require the imposition of boundary 
conditions along the outflow boundary T 0 (Streett & Macaraeg, 1989) zone technique. 

A similar treatment of the adjoint variables should have required consideration of 
an infinite domain [— oo<x<oo,0<y<h]. If this had been done, the boundary 
conditions (19) and (23) would not have been obtained along the inflow T^. In fact, the 
inflow boundary T^ for the state equation is the outflow boundary for the adjoint equations 
and, conversely, the outflow boundary T 0 for the state equation is the inflow boundary for 
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the adjoint equations. This is easily seen by comparing the leading inertial terms of the 
state and adjoint momentum equations (1), (2), (21), and (17), i.e. , with t increasing 


du 

~dt 


du du 

+ b v ~^~ 

dx dy 


and 


dv dv dv 

w + , ‘& +t % 


and with t decreasing 


du du du 

dt dx dy 


and 


dv dv dv 

dt dx dy 


Now, on both Ih and T 0 we have that u > 0 and v ~ 0 which is why is an inflow 
boundary and T 0 is an outflow boundary for the state. On the other hand, the fact that t 
is decreasing in the adjoint equations implies that now Ih is an outflow boundary and T 0 
is an inflow boundary for those equations. 

Thus, to be consistent with the treatment of the state equations, the adjoint outflow 
Ti should be treated in a manner similar to the above treatment of the state outflow T 0 . 
Thus, the boundary conditions (19) and (23) are not imposed along T^ and (21) and (17) 
are replaced by 


du „du dv du du dp 

dt dx dx dx dy dx 

, . d 2 u d 2 u . . _ ^ 

-vb(x)——-v—— = 0 m (0,T)xO c 


dx 2 


dy 2 


( 31 ) 


and 


dv ^du dv dv dv dp 

— — b u — — b^ — u w v — — 

dt dy dy dx dy dy 

d 2 v d 2 v 

-vb(x)—- - v—~ = 0 in (0,T)xfi c , 


dx 2 


dy 2 


( 32 ) 


respectively, where b(x) is a smooth function that is unity in the bulk of the flow, that 
vanishes in a neighborhood of the adjoint outflow boundary Tp and which is smooth 
throughout the flow. This treatment of the adjoint outflow does not require the imposition 
of any boundary conditions for the adjoint variables along r^. Finally, since T 0 is an inflow 
boundary for the adjoint equations, one should specify boundary conditions on the adjoint 
variables along that boundary segment. It can be seen that the proper conditions are given 
by 

u = 0 and 7 = 0 on (0,T) x T 0 . (33) 
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3.5. The Optimality System for Channel Flow 


We now have in hand the full optimality system for channel flow whose solutions 
determine the optimal states, controls, and adjoint states. These are given by (3)-(7), 
(14), (15), (18)-(20), and (22)-(33), where in (19) and (23) we do not impose the boundary 
conditions along IV Since (20) and (24) merely serve to determine the uninteresting 
Lagrange multipliers S 2 si, respectively, they can be ignored. Re-ordering and gathering 
the remaining equations, posed on the computational domain fie, yields the optimality 
system 


du du du dp . ,d 2 u d 2 u . ^ 

W- + M— + u— + - va(x)—-r-v—-r = 0 m (0, T) X fie 

at ox ay ox ox z dy z 


(34) 


dv dv dv dp d 2 v d 2 v 

a + “& + W + va{x) d. v - v df- = 0 


in (0, T) x fie 


du dv 
dx dy 


0 in (0, T) x fie , 


(35) 

(36) 


t = 0 


u 0 
^0 


in fie , 


(37) 



in (T 0 , Ti) 

in (0, T 0 ) and (Ti, T) , 


(38) 


u 


Ui 

Vi 


in (0,T), 


(39) 


m 


(0.T), 


du ^ du ^ dv du du dp d 2 u d 2 u 

~^ + u d^ + v d^~ u d^~ v d^~ d^~ u ^'dx 2 ~ V ~dy 2 

dv ^ du ^ dv dv dv dp d 2 v d 2 v 

~dt + “dy + V dy~ U di " - X) dV ~ ‘"By* 

|4 + |4 =0 in (OJIxD, 

dx dy 


t=T 


in fi 


c 


(40) 

0 in (0, T) X Sic , (41) 
0 in (0, T) X Sic , (42) 

(43) 

(44) 
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on 


' (o,r) x r 0 
(o, t) x r a 
< (o, t) x r w \r s 
(o, T a ) x r s 
K (Tb,T) x r s , 


a i i/ 


, du_ 

dy 


r i 


a 2 (~p + 2u^ - t 2 


on (T a ,T b ) x T a , 


(45) 


(46) 


d 2 gi 
dt 2 


+ 9 1 


1 / du 

JXdy 


on (T 0 ,Ti) x T a 


( 47 ) 


9 i |(=T 0 = 910(1) 


and 


dg 1 

at 


= 0 

t=Tx 


on r n 


(48) 


a ^2 

at 2 


+ 52 


1 

% 


p -\-2v 


dv 

dy 


on (T 0 ,Ti) x T a 


(49) 


and 


92\t=T 0 = 92 o(x) and 


dg 2 

at 


= 0 

t=Tx 


on 


r„ 


( 50 ) 


where r w denotes, in the channel flow case, the upper and lower boundaries of the finite 
channel except for the segment T a . 

The state equations (34)-(40) are driven by the given initial velocity (rto, r>o), the given 
inflow velocity (ui 7 Vi) } and the controls (gi } g 2 ). Indeed, the purpose of this study is to 
determine g\ and g 2 that optimally counteracts instabilities created upstream of T a . The 
adjoint equations (41)-(46) are homogeneous execept for the boundary condition (46) along 
T s , the part of the boundary along which we are trying to match the stresses. The data in 
that boundary condition is exactly the discrepancy between the desired stresses t\ and t 2 
and the stresses vdu/dy and — p + 2 vdu/dy along T s , weighted by the factors oq and a 2 . 
The equations for the controls (47)-(50) are driven by the negative of the adjoint stresses 
along T a , the part of the boundary along which we apply the control, weighted by the 
factors 1/fii and 1/ fi 2 . Of course this division into equations for the state, the adjoint 
state, and the control is really obscured by the fact that (34)- (50) are all intimitately 
coupled. 
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3.6. The Optimality System for Boundary-Layer Flow 

Following a similar process to that used in Sections 3. 1-3.5 for the channel flow case, 
one may derive an optimality system for the boundary layer flow case. The only difference 
is that in the latter case T^, denotes only the lower boundary with T a excluded and that 
the additional boundary condition (9) along the upper boundary T e must be accounted 
for. 

With the new interpretation for T w , one can still define the Lagrangian functional 
(12) and use the constraints (13) on allowable variations; however, due to (9), allowable 
variations are further constrained by 


Hr e 


Sp — 2v 


dSv 

dy 


r e 


0 for (0, T) . 


(51) 


Note that (51) implies that, along T e , one may not choose the variations in Sp and dSv/dy 
independently. Considering, simultaneously, variations in p, v, and dv / dy along T e , one 
obtains 


»T 


f ( dSv 

/ v ( Sp — 2v — — ) dYdt 

<r e V dy 

>T 


+ / / Sv fp + 2z/ — + vv ) dYdt 

o r g \. y 


T 


„ dSv 

u — — aT dt 
Ox 


0 . 


(52) 


The first integral in (52) vanishes due to (51). One can show that 


u = 0 on (0,T) x T e (53) 

so that the third integral in (52) vanishes as well. Then, letting Sv be arbitrary along T e 
in (52) yields 

p + 2^^— + vv = 0 on (0, T) x T e . (54) 

The resulting optimality system for the boundary layer flow case is given by (9), 
(34)-(50), (53), and (54), where in (40) and (45) T^ notes the finite computational lower 
boundary with T a excluded and in (53) and (54) T e denotes the finite computational 
uppper boundary. 
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4. NUMERICAL EXPERIMENTS 


Here, the optimal control methodology developed in Section 3 is applied to a boundary 
layer flow having a single instability wave that can be characterized by a discrete frequency 
within the spectrum. We are not concerned with the details of how disturbances are 
ingested into the boundary layer; the underlying assumption here is that natural transition 
involves some dominant disturbances that can be characterized by waves, and in fact, 
in the present study, by a single wave. (In a subsequent study, controlling transitions 
which consist of unsteady, three-dimensional instabilities will be explored.) As described 
in Joslin et al. (1995) (among others), these discrete small-amplitude instabilities can be 
suppressed through wave cancellation (WC) using known exact information concerning the 
wave. Hence, the optimal control is “known” for validation of the present DNS/optimal 
control theory numerical approach in which the instability is to be suppressed without any 
a priori knowledge of said instability. 

The formidable coupled system (9), (34)-(50), (53), and (54) is solved in an iterative 
manner. First, a guess is made for the controls gq and g 2 (typically one starts with no 
control, i.e., gi = g2 = 0) and then the sub-system (9) and (34)-(40) is solved for the 
state variables, i.e., the velocity field (u, v) and pressure p. Then, using these velocity and 
pressure fields, the sub-system (41)-(46), (53), and (54) is solved for the adjoint or co-state 
variables (u,0) and p. Then, using these adjoint variables, the controls gq and <72 are then 
found by solving the sub-system (47 )- (50) . The procedure is repeated until satisfactory 
convergence is achieved. 

The nonlinear, unsteady Navier-Stokes equations and linear adjoint Navier-Stokes 
equations are solved by direct numerical simulation (DNS) of disturbances that evolve 
spatially within the boundary layer. The spatial DNS approach involves spectral and 
high-order finite-difference methods (Joslin et al 1992, 1993) and a three-stage Runge- 
Kutta method (Williamson 1980) for time advancement. The influence-matrix technique 
is employed to solve the resulting pressure equation (Danabasoglu, Biringen & Streett 1991 
and Street & Hussaini 1991). Disturbances are forced into the boundary layer by unsteady 
suction and blowing through a slot in the wall. At the outflow boundary, the buffer-domain 
technique of Streett & Macaraeg (1989) is used. 

In the present study only normal injection or suction control is allowed, so that we 
set g\ = 0 in (38), = 0 in the functional (11), and ignore (47) and (48). Also, we only 

match the normal stress along T s so that we choose oq = 0 in the functional (11) and in 

(46), 
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4.1. Computational Parameters 

For the computations, the grid has 401 streamwise and 41 wall-normal points. The 
free stream boundary is located 756* from the wall, and the streamwise length is 2246* 
which is equal to approximately 8 TS wavelengths. The nondimensional frequency for 
the forced disturbance is F = u / R X 10 6 = 86; the forcing amplitude is Vf = 0.1%. 
The Reynolds number based on the inflow displacement thickness (6*) is R = 900. (The 
boundary segment along which disturbance forcing and control is effected as well as where 
stress matching occurs are located within the unstable region of the linear stability neutral 
curve.) A time-step size corresponding to 320 steps per period T p is chosen for a three- 
stage Runge-Kutta method. Based on the disturbance frequency, a characteristic period 
can be defined as T p = 2tt / u = 81.1781; the resulting time-step size is then At = 0.2537. 

To complete one period of the active-control simulation process, 0.75 minutes on the 
Cray C-90 are required using a single processor. Note, two periods of cost (T a — > T& and 
T& — > To) are required to complete one iteration of the DNS/adjoint system. Although in 
general any time interval may be specified for T a — > Tj, this study uses integer increments 
of the period (T p ) for simplicity. Hence, T a — > T& = 2 T p would cost 4 T p in computations, 
or roughly 3 min of C-90 time per iteration. Because only a single small-amplitude wave 
(linear) is forced, the above grid is more than adequate; however, a grid refinement was 
performed and produced results equivalent to the results reported here. 

For this study, the disturbance forcing slot Tj, the control or actuator orifice T a , 
and the matching or sensor segment have equal length 4.48 6*. The forcing is centered 
downstream at 389 . 62 A* (the Reynolds number based on the displacement thickness at 
that location is R = 1018.99), the actuator is centered at 403.62 6* (R = 1037.13), and the 
sensor is centered at 417.62 6* ( R = 1054.97). These separation distances were arbitrarily 
chosen for this demonstration. In practice, the control and matching segments should have 
a minimal separation distance so that the pair can be packaged as a single unit, or bundle, 
for distributed application of many bundles. 

4.2. Results 

All simulations allow the flow field to develop for one period, i.e., from t = 0 — > T a = T p 
before control is initiated. In the first series of simulations, the interval during which control 
is applied is arbitrarily chosen to be T a — > T& = 2 T p . Based on oq = /% = 0, «2 = 1, and 
@2 = 10, the wall-normal velocity and disturbance energy results are shown in Figure 2 
for the series of iterations leading to convergence. Figure 3 shows the convergence history 
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for the actuator g 2 and measured normal stress — p + 2 vdv / dy as a function of time. The 
velocities are obtained at a fixed distance from the wall corresponding to 1.186* and the 
measure of energy is given by E = u 2 + v 2 ; both velocity and energy are obtained at the 
fixed time T 5 . Convergence is obtained with 8 iterations; however, for all practical purposes, 
the results with 4 iterations could be used. The results demonstrate that a measure of 
wave cancellation can be obtained from the DNS/control theory system. The wall-normal 
amplitude of the modified wave at R = 1092.5 is 40 percent of the uncontrolled wave; the 
control without optimizing the choice of op, « 2 , /?i , and f3 2 has led to a 60 percent decrease 
in the amplitude of the travelling wave. Clearly, Figure 2 shows that a net reduction of 
the disturbance energy is obtained by energy input due to the control. This results in a 
delay of transition by-way-of a suppression of the instability evolution. 

In the simulation, the control has been applied from T a — > T& only; therefore, for 
t > T 5 , (5) or (38) indicate that the actuation is discontinued. Figures 4 and 5 compare 
the converged results (Cl= 8 th iteration of Figures 2 and 3) with results for one period 
after control t = T& + T p . The measured disturbance tends toward the uncontrolled so- 
lution when the actuation is discontinued (as expected); because the control was applied 
for T a — T p — > T}, = 2 T p , 2 periods are required after T& before the computed solution 
in the window exactly matches the uncontrolled solution. The profiles of uncontrolled, 
control, and discontinued control flows are shown in Figure 5 at the downstream location 
corresponding to R = 1073.2. Clearly, the control only removes energy from the sys- 
tem (decreases the wave instability amplitude); the resulting profiles retain the expected 
instability profile shape. 

The effect of varying the window size (T a ,T&) is shown in Figures 6 and 7. The 
previous converged results (Cl) are shown with converged results (C 2 ) for the extended 
window (T a = T p — > T& = 3 T p ). The results are identical for the first two periods of 
time and indicate that extending the amount of time for control serves to extend control 
only. This result also indicates that one can solve for the optimal control over a given time 
interval (T a ,T&) by breadking up that inerval and solving for the optimal control over a 
series of smaller subintervals. This approach usually leads to substantial savings in CPU 
and memory costs. 

Figure 7 reveals an additional insight about the present DNS/control theory. The re- 
sulting optimal control g 2 approaches the desired wave-cancellation time-periodic solution 
as the temporal length (T a ,T&) is increased. This is convincing evidence that the present 
self-contained methodology is valid. 
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Figures 8 and 9 show velocity contours of the uncontrolled and controlled (C2) insta- 
bility waves. Both results indicate that the instability wave is growing with downstream 
distance (left-to-right); however, the C2 contours are significantly reduced in magnitude. 
The furthest downstream levels of C2 approach the level of the uncontrolled wave at the 
furthest upstream location. If the growth rates are assumed to be the same, then a three- 
wavelength transition delay has been achieved by applying the optimal control. 

The instability wave resulting from wave-cancellation (WC) is shown with the control 
(C2) in Figures 10 and 11. For the present comparison, the amplitude of the actuation 
for WC was adjusted until nearly exact wave cancellation was achieved. Although the 
DNS/control theory did not achieve the same level of energy removal, the similar effect of 
WC was achieved without any a priori knowledge of the instability. Also, note that Figure 
11 shows that the optimal control of the control theory has nearly the exact phase charac- 
teristics as WC and only lacks the necessary amplitude for additional wave cancellation. 
These encouraging results suggest that by the appropriate selection of op, « 2 , /?i , and f3 2 , 
the optimal control can be made nearly as effective method of instability suppression as 
exact wave cancellation. 

From the wave-cancellation study of Joslin et al (1995), the relationship between 
amplitude of the actuator (v a ) with resulting instabilty can be shown in Figure 12. A 
similar result was shown in the channel flow wave-cancellation study in Biringen (1984). 
The trend indicates, that beginning with a small actuation amplitude, as the actuation 
level is increased, the amount of wave cancellation (energy extraction from the disturbance) 
increases. At some optimal actuation, nearly exact wave cancellation is achieved for the 
instability wave. As the actuation amplitude further increases the resulting instability 
amplitude increases; this was clearly explained in Joslin et al (1995) to occur because in 
the wave superposition process, the actuator wave becomes dominant over the forced wave. 
At this point, the resulting instability undergoes a phase shift corresponding to the phase of 
the wave generated by the actuator. The relationship depicted in Figure 12 is encouraging 
for the DNS/optimal control theory approach and suggests that a gradient descent type 
algorithm might further enhance the wave suppression capability of the present approach. 
Namely, an approach for the optimal selection of op, « 2 , /?i , and @2 might lead to a more 
useful theoretical/computational tool for flow control. 

To simply demonstrate this concept, Lagrange interpolation (or perhaps extrapola- 
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tion) is introduced for f3 1 and f3 2 based on imposed values for a 1 and 0 / 2 : 


/3d 1 = 


PlM,2 T l,2 ) ft 1,2 ( T l*,2 r l) 2 ) 

( T U2 — r l,2 ) 


(55) 


where r* 2 are some desired values of the stress components and t™ 2 are the stress com- 
ponents based on the choice /?” 2 . Although 7 % and r 2 * may be equivalent to the target 
values T\ and T 2 in the functional (11), this may lead to significant over/under shoots for 
the interation process. Instead, r * and r 2 * is the incremental decrease, or target value, for 
interpolation to more desirable /% and f3 2 values. To illustrate this process, the f3 2 = 10 
(C2) and /?2 = 11 control results are obtained with the iteration procedure. The measures 
of normal stress are somewhat arbitrarly obtained at some time as measured by the sensor 
or matching segment T s ; the values of the normal stress are given in the Table 1. These 
values are used for a desired normal stress r 2 , which in this case is 65% of the @2 = 11 
results. 


h 

normal stress 

10 

11 

9.369 x 10“ 6 
8.814 x 10“ 6 


Table 1. Normal stress for two values of /%. 


Using the results for f3 2 = 10 and f3 2 = 11 in (55) yields the value f3 2 = 16.5 which 
is used in a simulation to obtain a greater degree of instability suppression. The WC 
results and the enhanced optimal control (C3) solution are shown in Figures 13 and 14. 
This interpolation approach based on relationship of Figure 12 indicates that optimizing 
f3 2 has led to results very close to WC. The solutions differ somewhat near t = T a and 
t = Tb because of the conditions (48) and (50) that serve to control the levels of gq and 
g 2 - For all practical purposes, the solutions obtained with the present DNS/control theory 
methodology yield the desired flow control features without prior knowledge of the forced 
instability. 

The adjoint system requires that the velocity field (u } v) obtained from the Navier- 
Stokes equations (34)-(40) be known for all time. For the iteration sequence and a modestly 
course grid, 82 Mbytes of disk (or runtime) space are required to store the velocities at 
all time steps and for all grid points. For T a — > Tj, = 3 Tj,, 246 Mbytes are necessary 
for the computation. Clearly for three-dimensional problems the control scheme becomes 
prohibitively expensive. Therefore, a secondary goal of this study is to determine if this 
limitation can be elliminated. 


25 



Because the characteristics of the actuator (gq and gq) and resulting solutions are 
comparable to WC, some focus should be placed on elliminating the enormous memory 
requirements discussed above. This limitation can easily be removed if the flow-control 
problem involves small-amplitude unsteadiness (or instabilities). The time-dependent coef- 
ficients of the adjoint system (41)-(42) reduce to the steady-state solution and no addition 
memory is required over the Navier-Stokes system in terms of coefficients. This has been 
verified by a comparison of a simulation with steady coefficients compared with the C2 
control case. The results shown in Figure 15 are identical (as expected). Additionally, 
if the instabilities have small amplitudes, then a linear Navier-Stokes solver can be used 
instead of the full nonlinear solver, which was used in the present study. This linear system 
would be very useful for the design of flow-control systems. However, if the instabilities in 
the flow have sufficient amplitude to interact nonlinearly, then some measure of unsteady 
coefficient behavier is likely required. Depending on the amplitudes, the coefficients saved 
at every time-step may be replaced with storing coefficients every 10 or more time-steps 
thereby reducing the memory requirements by an order of magnitude. This hypothesis 
will require validation in a future study. 

CONCLUSIONS 

The coupled Navier-Stokes equations, adjoint Navier-Stokes, and optimality condition 
equations were solved and validated for the flow-control problem of instability wave sup- 
pression in a two-dimensional, flat plate, boundary layer. By solving the above system, 
optimal controls were determined that met the objective of minimizing the perturbation 
normal stress along a portion of the bounding wall. As a result, the optimal control 
was found to be an effective means for suppressing two-dimensional, unstable Tollmien- 
Schlichting travelling waves. The results indicate that the DNS/control theory solution 
is comparable to the wave-cancellation result but, unlike the latter, requires no a priori 
knowledge of the instability characteristics. 
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A - Actuator 
S - Sensor 


C 


C - Optimal Control Theory 


Figure 1. Schematic of active flow control using optimal control theory. 
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Figure 5. Disturbance velocity profiles for no control ( T =4 ) , contro 1(C1), and after control 
is used and turned-off (C=4) in flat-plate boundary-layer flow. 
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Figure 6. Disturbance velocity with downstream distance for no control ( T =4 ) and control 
(Cl for Ti — To = 2 T p and C2 for T\ — To = 3 T p ) in flat-plate boundary-layer flow. (Velocity 
signal at y = 1.18(5* from wall.) 








Figure 8. Contours of (a) u and (b) v velocities for no control in flat-plate boundary-layer 
flow. (Contours: 7.5 X 10 -5 to —7.5 X 10 -5 ; intensity increases left-to-right ) 
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(b) 


Figure 9. Contours of (a) u and (b) v velocities for control in flat-plate boundary-layer 
flow. (Contours: 7.5 x 10 -5 to —7.5 x 1 0 — 5 . ) 
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Disturbance velocity and energy with downstream clist; 
rol (C2), and wave cancellation (WC) in flat-plate bound 
at it = 1.18(5* from wall; Ti — Tn = 3 T n .) 









Figure 12. Disturbance velocity resulting from variations in actuator amplitude from 
simulations in Joslin it al (1995). 













